Stable Distributions in Stochastic Fragmentation 
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We investigate a class of stochastic fragmentation processes involving stable and unstable frag- 
ments. We solve analytically for the fragment length density and find that a generic algebraic 
divergence characterizes its small-size tail. Furthermore, the entire range of acceptable values of de- 
cay exponent consistent with the length conservation can be realized. We show that the stochastic 
fragmentation process is non-self-averaging as moments exhibit significant sample-to-sample fluctu- 
ations. Additionally, we find that the distributions of the moments and of extremal characteristics 
possess an infinite set of progressively weaker singularities. 
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I. INTRODUCTION 

Ftagmentation is a basic stochastic process with a va- 
riety of applications ranging from geology ^ and frac- 
ture 1^ to the breakup of liquid droplets and atomic 
nuclei and impact fragmentation of solid objects 

1^-^. Fragmentation processes are relevant to numer- 
ous seemingly unrelated problems such as spin glasses, 
Boolean networks and genetic populations [p|-p^. This 
study is particularly motivated by computational biol- 
ogy applications involving partitioning of nonstationary 
time series into stationary sub-series, and conversely, re- 
construction of a series from its sub-series. Examples 
include "parking" strategies for genome reconstruction 
1^-0 and DNA segmentation algorithms 

The primary feature of the fragmentation process 
which underlies segmentation algorithms is that rather 
than continuing indefinitely, the process stops when a 
segment reaches an acceptable "homogeneity" level. Fur- 
thermore, fragmentation does not continue indefinitely in 
impact fragmentation processes where the kinetic energy 
and the size of the fragments determine the number of 
fragmentation events. These applications lead to stochas- 
tic fragmentation processes where not only the way by 
which fragments are produced, but also the number of 
fragmentation events is subject to a random process. 

In the simplest stochastic fragmentation process, one 
starts with a fragment and breaks it into two pieces. 
With probability p, a newly formed fragment remains 
unstable, i.e, it continues to participate in fragmenta- 
tion events, while with probability q = 1 — p, it be- 
comes stable, and is never fragmented again. The pro- 
cess is repeated for all unstable fragments until all frag- 
ments become stable. For this simple process, we show 
that the final length density, P{x), is purely algebraic, 
namely P{x) = 2qx~^^ [Q. Similar scale free behavior 
was observed in other fragmentation processes |p2| , p3| . 
Stochastic fragmentation also exhibits intriguing statis- 
tical characteristics including moments which are non- 
self-averaging, essential singularities in the distribution 
of these moments, and an infinite set of singularities in 



the distribution of the largest and the smallest fragments. 

In this paper, we study statistical properties of stochas- 
tic fragmentation processes. In Sec. H, we introduce the 
stochastic fragmentation model and derive the fragment 
length density P{x). We also consider generalizations 
to size dependent fragmentation densities and fragmen- 
tation probabilities, and obtain the exponent underlying 
the algebraic behavior of P(x) as a root of a transcen- 
dental equation. In Sec. Ill, we show that the stochastic 
fragmentation process is non-self-averaging. Specifically, 
the moments Yq = xf exhibit significant sample-to- 
sample fluctuations. In Sec. [V, we study extremal char- 
acteristics such as the distribution of the largest and the 
shortest fragment, and we find that both the extremal 
distributions and the distribution of the moments are 
characterized by an infinite set of progressively weaker 
singularities. In Sec. we present an application to ran- 
dom sequential adsorption processes, and we summarize 
this work in Sec. 



II. THE MODEL 

In the basic stochastic fragmentation process, we start 
with the unit interval, which is considered to be unsta- 
ble. This interval is fragmented into two pieces of length 
I and 1 — ^, where I is drawn from a uniform probability 
density p{l) — 1. With probability p, each of these frag- 
ments remains unstable, while with probability g = 1 — p, 
it becomes stable and does not undergo further fragmen- 
tation. The process is iterated for unstable fragments 
until all fragments become stable. 

The average total number of stable fragments, (iV), 
can be directly evaluated. Consider a fragment pro- 
duced in the first fragmentation event. With probability 
q it is stable, and consequently, only a single fragment 
is produced; otherwise, the process is repeated. Hence 
(TV) = 2((z + p(7V)), yielding 



(iV) 



r2g/(l-2p), 
1 oo, 



ifp< 1/2; 
if p> 1/2. 



(1) 
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The average total number of fragments diverges as the 
probabiUty p approaches the critical point pc = 1/2, re- 
flecting the critical nature of the corresponding branching 
process j24j. 

Next, consider the fragment length density, P{x), of 
stable fragments of length x. The recursive nature of the 
process can be used to obtain the governing equation for 
the fragment length density 
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(2) 



The first term accounts for stable fragments produced at 
the first generation, and the second term describes cre- 
ation of an .T-fragment from a larger y- fragment. 

Equation can be solved by employing the Mellin 
transform technique. Let 



M{s) = / dxx^'^Pix). 



(3) 



Hereinafter, integration with unspecified limits is car- 
ried over the unit interval. Equations (^) and yield 
M{s) = 2s^^ [q + pM{s)\, and consequently, the Mellin 
transform of the length density reads 



M{s) 



2q 



2p 



(4) 



Note that the total length is conserved, M(2) = 1, and 
that the total number of stable fragments, M(l) = (N), 
is consistent with Eq. (|l|). 

The length density can be obtained by inverting the 
Mellin transform 



P{x) ^2qx 



-2p 



(5) 



Remarkably, the length density is purely algebraic over 
the entire range < a; < 1. Generally, given an al- 
gebraic divergence near the origin, P{x) ^ x '', length 
conservation provides the upper bound 7 < 2, and as 
< 7 = 2p < 2, the entire range of possible divergences 
is realized by tuning p. 

Interestingly, scale free distributions were also found 
for a dual stochastic aggregation process where aggre- 
gates may turn stable after each aggregation event. In 
that case, the large size tail of the distribution decays al- 
gebraically 1^ . We also note that algebraic distributions 
have been observed in a number of recent impact frag- 
mentation experiments involving rods, spheres, bricks, 
etc., with the corresponding decay exponents typically 
ranging between 1 and 2 

Notice that the form of the fragment length density is 
not altered as the critical point pc = ^ is passed. Never- 
theless, this point is characterized by a unique property. 
Starting from an interval of length Lqj Eq. (^) can be 
generalized to yield 



Pix) = 2q^ 



2p- 



(6) 



Hence, at the critical point, P{x) becomes independent 
of the initial interval length Lq. 

In the above basic model both the fragmentation den- 
sity and the fragmentation probability were independent 
of the fragment length. In the following, we show that 
even when these functions become length-dependent, 
P{x) remains algebraic in the small size limit. 



A. Arbitrary fragmentation density p{l) 

Consider a fragmentation process in which an interval 
is broken into two fragments of relative lengths I and 1 — I 
with an arbitrary fragmentation density p{l). This den- 
sity satisfies the normalization constraint J dl p{l) — 1 
and the symmetry requirement p{l) = p{\ — I). The gov- 
erning equation for the fragment length density reads 



P{x)^2qp{x) + 2p I ^p(y)p(^ 

V \V 



(7) 



The Mellin transform (0) of the length density satisfies 
M{s) — 2p{s) [q +pM(s)\, and consequently 



M{s) 



2q 



^i-^{s)-2p' 



(8) 



where p{s) = J dll''~^p{l) denotes the Mellin transform 
of the fragmentation density p{l). The symmetry of the 
fragmentation density implies /i(2) = 1/2 and therefore 
M(2) = 1, which confirms the conservation of length. 
Also, the normalization condition implies /i(l) = 1, 
which confirms that the average total number of frag- 
ments is given by Af (1) = (TV), in agreement with Eq. (]^). 

The Mellin transform shows that the fragment length 
density is scale free only when the fragmentation density 
is uniform. Nevertheless, the algebraic small-size behav- 
ior remains robust. Indeed, Eq. (||) suggests that the 
most important property of M{s) is a simple pole whose 
location s = 7 is found from relation 2pp(^) = 1. This 
simple pole implies a power-law asymptotics of the frag- 
ment length density 



P{x) ~ Ax-T 



(9) 



as a; — > 0. The exponent 7 can be determined from the 
entire fragmentation density p{x) via the relation 



2p I dir'^p{l) = 1. 



(10) 



The prefactor in Eq. (||) reads A = [qp{'^)\/[pp,' {'^)], 
which is simply the residue of the pole at 5 = 7. 

At the critical point Pc = 1/2 one has 7 = 1, indepen- 
dent of the fragmentation density p{l). The relation ( p^ ) 
also shows that a generic behavior 7^2 occurs if the 
probability of becoming stable vanishes, i.e., if p —^ 1. 
In the complementary p — > limit, the small-size behav- 
ior of p{l) determines the small fragment distribution. In 



2 



particular, if p{l) in the limit I 0, Eq. ( [ic| ) shows 

that 7 ^ r as p ^ 0. Hence, in this case the restricted 
exponent range r < 7 < 2 emerges by tuning p. 

As an illustration, consider the fragmentation density 
p{l) = B[1{1 - l)]^-^, with B = T{2S)/r^{S) ensuring 
proper normalization. The exponent 7 is determined 
from Eq. ( |lO| ) to give 



2p 



r(2(5)r(7 + J - 1) 

T{5)T{-f + 2d-l) 



= 1. 



(11) 



This relation shows that the exponent 7 always belongs 
to the range 1 — (5 < 7 < 2. In the extreme case of ^ ^ 0, 
the decay exponent is concentrated near 7=1. Such 
universal P{x) ~ behavior is empirically observed in 
DNA segmentation algorithms . In the other extreme 
(5 — *■ 00, the exponent simplifies to 7 = 1 + ln2p. Note 
also that explicit results for both P{x) and the exponent 
7 can be obtained for integer ^'s. The case of (5 = 1 cor- 
responds to the uniform density. For 6 = 2, the Mellin 
transform reads ^(s) = 6/(s-|-l)(s-l-2), and Eq. yields 
the following fragment density 



Fix) 



Uq 



3-Vl + 48p 



X ^ J.. (12) 



Generally, the length density is a linear combination of 6 
power laws for all integer 6's. 



B. Arbitrary fragmentation probability p{x) 

We now discuss the complementary generalization, in 
which the probability p{x) that a new fragment remains 
unstable depends on the fragment size x. This may ac- 
tually be the case in impact fragmentation as well as 
in DNA segmentation, where fragments have an intrin- 
sic size scale below which the fragmentation probability 
becomes negligible. For an arbitrary fragmentation prob- 
ability p{x) the governing equation reads 



P{x) = 2 



l-p{x)+ r^p{y)p(^ 
Jx V \V 



(13) 



Consequently, the Mellin transform of the length density 
admits the general solution 



^ ' si- 2a{s) 



(14) 



where (t[s) = J dx x^^^p{x) is the Mellin transform of 
the probability p{x). 

The small size tail is determined by the poles of M(s). 
When the condition 2a{0) > 1 is satisfied, M{s) has a 
simple pole at s = 7, and therefore, the small size behav- 
ior remains algebraic as in Eq. (^. The corresponding 
exponent 7 is determined from 2cr(7) = 1, or explicitly 



dx x"' ^p{x) = 1, 



(15) 



and the prefactor A = Xr2]/[27cr'(7)] equals the residue 
at s = 7. Equation (15|) is a transcendental equation, 
and the details of the function p{x) in the entire range 
< X < 1 determine the exponent 7. This situation is 
reminiscent of the behavior found when the fragmenta- 
tion probability along the interval was not uniform. 

In the complementary case of 2a{0) < 1, Eq. (Ill) 
shows that M{s) has a simple pole at s = 0, viz. M{sj = 
2s"V[1-2(t(0)] + . . . This implies that the length density 
is regular in the small size limit, P{x) — > 2/[l — 2cr(0)] 
as a; — > 0. In the marginal case 2a{0) = 1, we find 
M{s) ~ s~^, which leads to a logarithmic divergence of 
the length density, i.e. 

As an illustration, 
p{x) = x^ with A > 0. In this case a{s) = l/(s -I- A), 
and consequently M(s) = 2A/[s(s-|-A — 2)]. Inverting the 
Mellin transform M{s) gives the fragment length density 



P{x) ~ In J as a; ^ 0. 

consider the solvable example 



P(x) 



f a: - 1) 

41nl- 

X 

X-2 



< A < 2; 
A = 2; 
2 < A. 



(16) 



Again, the range < 7 < 2 becomes accessible. 

Hence, the algebraic small size divergence is robust as 
it extends to situations where either the fragmentation 
density or the fragmentation probability are size depen- 
dent. The entire form of these functions is needed to 
calculate the corresponding power-law exponent. 

In the rest of this paper, we restrict ourselves to the 
basic model where the fragmentation density is uniform, 
p{x) = 1, and the fragmentation probability p{x) = p is 
size independent. 



III. THE MOMENTS 



The fragment size distribution represents an average 
over infinitely many realizations of the stochastic frag- 
mentation processes, and hence, it does not characterize 
sample-to-sample fluctuations. In this section, we show 
that fluctuations do not vanish in the thermodynamic 
limit, and therefore, the process is non-self-averaging. We 
investigate sample-to-sample fluctuations by computing 
the moments Yq defined by 



(17) 



where the sum runs over all fragments in a given realiza- 
tion. These moments have proved useful in a variety of 
contexts including spin glasses, random maps, and ran- 
dom walks ||2^,|2^. The fact that these moments have 
non-trivial probability distributions is a signature of lack 
of self-averaging. 

Before we attempt to derive these probability distri- 
butions, we start with the simpler task of computing the 
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expected values of the moments (Fq) as well as their cor- 
relations (XoXis)- For integer a, (Yq) is the probabihty 
that a randomly chosen points in the unit interval be- 
long to the same fragment. Similarly, for integer a and 
integer /3, (YaYa) is the probability that among a + (3 
points chosen at random, the first a points all lie on the 
same fragment, and the last [3 points all lie on another 
(possibly the same) fragment. 
The expected value of satisfies 



(18) 



The g-term corresponds to the situation where a first gen- 
eration fragment becomes stable while the second term 
describes the complementary situation. Equation ( p^ ) 



gives 



(Y^ 



2q 



l-2p 



(19) 



if a > 2p—l, and (Fq) = cx) if a < 2p—l. Since the single 
point averages are simply the moments, (Yq) = M(a-\-l), 
Eq. ( p^ ) indeed agrees with Eq. (||). 

Higher-order averages cannot be computed from the 
fragment size density. However, one can obtain exact 
expressions for higher averages from recursion relations 
similar in spirit to Eq. (p^). For instance, (YaFa) satisfies 

{Y^Yf}) = 2{q + p{Y^Y0)) J dyy''+^ (20) 

+ 2iq + p(Y^)){q + p{Yp)) f dyy'-il-yf. 



(21) 



The solution to the above recursion relation reads 



{Yc^Yp) 



C[a,P){q+p{Y^)){q+p{Yp)) 
a + P + l-2p 



if a, (3,a + (3 > 2p — 1, and (YaYg) — oo otherwise. Here 
we used the shorthand notation C(a,/3) = ^^Y{a+p+i)^^ ■ 

Equation (|l|) shows that {YaYf}) ^ (Fa)(Ys) and, 
in particular, (Y^) ^ (YaY ■ Hence, fiuctuations in Yq, 
do not vanish in the thermodynamic limit, which im- 
plies that the stochastic fragmentation process is non- 
self-averaging. This means that statistical properties ob- 
tained by averaging over all realizations are insufficient 
to probe sample-to-sample fluctuations. Lack of self- 
averaging was also found in fragmentation processes that 
exhibit a shattering transition p9|-|3l|| . 

It is possible to evaluate higher-order averages such as 
(Y^). However, even for small n these averages become 
quite cumbersome and not terribly illuminating. Instead, 
one might try to obtain the distribution, QaiYa), of pos- 
sible outcomes of the moments Yq. Let us first consider 
the fragment number distribution Qo{N) (the zeroth mo- 
ment equals the number of fragments, Yq = A'^), which 
can be determined analytically. The minimal number of 
fragments is produced when both of the first generation 



fragments are stable, and hence, Qo{2) = q^ . Similarly 
for > 3 we obtain the recursion relation 



Qo{N) = 2pqQoiN 

Ni+N2=N 



1) (22) 
Qo{Ni)Qo{N2), 



where the total number of fragments N is obtained in 
various ways from a smaller number of fragments that 
appear after fragmentation of the two first generation 
fragments. Specifically, if exactly one of the first genera- 
tion fragments is unstable, it should produce A^ — 1 stable 
fragments. If both of the first generation fragments are 
unstable, they can produce Ni and N2 fragments, re- 
spectively, subject to the constraint A^i + iV2 = N. This 
explains the right-hand side of Eq. (^) . 

Equation ( p2| ) can be solved by introducing the gen- 
erating function Qq{z) = J2N>2Qoi^)^^ ^ which satis- 
fies p'^Qliz) -t- (1 - 2pqz)Qo{z) + q'^z^ = 0. Solving this 
quadratic equation yields the generating function 



1 — 2pqz — ^1 — Apqz 



Expanding Qo{z) in powers of z gives 



Qo{N) 



r(i)r(iv + i) V 



(23) 



(24) 



At the critical point, pc = 1/2, the number distribution 
decays algebraically in the large A^ limit: 



(25) 



In the vicinity of pc = 1/2, the number distribution at- 
tains the scaling form 



Qo{N) - Ar-3/2g^p [_47v(Ap)2] 



(26) 



where Ap ^ pc — p. Hence, below the critical point, the 
tail of the number distribution is exponential. 

The probability that an infinite number of fragments 
is produced is given by 

oc 

Qo(oo) - 1- Qo(A^) = l-go(2 = 1). (27) 

N=2 

Below the critical point, the number of fragments re- 
mains finite, i.e., Qo{oo) = 0. In the complementary 
case of p > Pc, with finite probability, an infinite number 
of fragments is produced 



Qo(oo) = 1 2 ■ 



(28) 



The case a = is unique in the sense that the vari- 
able Yq = A^ is discrete. Another special case is a = 1, 
where length conservation dictates Yi = 1, and therefore 
the distribution is trivial, Qi{Yi) = S{Yi — 1). When 
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< a < 1, then > 1, and Qa has support in the in- 
terval (1, oo), as in the case of a = 0. Equation ( ^ ) then 
suggests that for p > Pc the distribution Qa{Ya) should 
have a singular component at Yq = c». 

In the following, we focus on the more interesting case 
of a > 1. Here, the inequality Yq < Yi = 1 implies 
that the distribution (5q(Yq) has support on the interval 
(0, 1). If both of the first generation fragments happen 
to be stable, then Y^ = a;" + (1 — x)", where x is cho- 
sen uniformly in the unit interval. Both fragments are 
stable with probability q^, and the corresponding con- 
tribution to Qa(Ya), which we denote by Ila{Ya), reads 
Hq {Ya) = 2q^ , where x is the greater of the two roots 
of equation Ya = x" + {1 — .x)" (the factor 2 accounts 
for the smaller root). Although it is generally impossible 
to express the above formula solely in terms of Ya, in 
some special cases one can determine Ila{Ya) explicitly; 
for example, 02(^2) = g^(2Y'2 - l)-!/^. 

Note that the distribution 112 (Y2) has a singularity at 
Y2 = 1/2, which obviously implies a singularity of <52(^2) 
at the same point. To understand the origin of this singu- 
larity, notice that when the process ends with two stable 
fragments, then I2 = + (1 — x)^ > 1/2. Therefore, the 
behavior of (52(5^2) for the case of Y2 < 1/2 is noi affected 
by realizations with two final fragments, and this explains 
the singularity at ^2 — 1/2. If the process ends with three 
stable fragments, then Y2 — x\+X2 + {\^ — xi — X2)'^ > 1/3. 
Similarly, if we end with k stable fragments, then Y2 > 
1/k. Hence, we anticipate that the distribution Q2{Y2) 
has singularities at I2 = 1/^ for integer k > 2. Similar 
singularities underlie distributions of moments in a num- 
ber of random processes, including random walks, spin 
glasses, random maps, and random trees 0,|ll|,||,|7|. 

A straightforward generalization of the above argu- 
ment suggests that for arbitrary a > 1, the distribution 
Qa{Ya) possesses singularities at Ya = k^^°'. The ex- 
istence of these infinitely many singularities shows that 
analytical determination of the distribution Qa{Ya) is 
hardly possible. Indeed, Qa{Ya) satisfies the difficult in- 
tegral equation 



.iYa)^IlaiYa)+2pq 



dl 



(1 



dZ 



dl 



^"(1 - 1)°" 



Qo 



- 0" 
z 



Qo 



Qc 



(1-0" 

Ya-Z 

(1-0" 



(29) 



Equation ( p9| ) has been derived by repeating the steps 
used in the derivation of Eq. (p2[). The first (second) 
term on the right-hand side of Eq. ( ^9|) corresponds to 
the case where two (one) of the first generation fragments 
are stable. The third convolution term describes the al- 
ternative case when both of the first generation fragments 
are unstable. Note that in addition to the recursive na- 
ture of the process, we have employed extensivity, i.e., 
(Ya) (x I", in an interval of length I. 

In order to study the small- Y^ behavior of the distri- 
bution, we employ the Laplace transform method. From 
Eq. (H), RaiX) = J^dYae-^^-QaiYa) obeys 



(30) 



In Eq. (pO|), we do not write explicitly the Laplace trans- 
form of the first two terms of Eq. (^9|) , because these two 
terms become negligible in case of Ya 0. The Y^ ^ 
asymptotics of Qa{Ya) is reflected by the A — > 00 asymp- 
totics of i?Q,(A). We argue that Ra{X) ~ exp(— ^A") 
with uj = 1/a. Assuming that aoj > 1, the above 
stretched exponential form of Ra{X) shows that the prod- 
uct Ra [XI"] Ra [A(l — 0"] would reach a maximum that 
greatly exceeds Ra{X) at I — 1/2, in contradiction with 
Eq. (pOj). Alternatively, if aoj < 1, the above prod- 
uct would reach its maximum at ^ = and 1 = 1. 
Then, the integral on the right-hand side of Eq. ( |30| ) 
would become 2Ra{\) fdl Ra [M"] oc \-^/°'Ra{X), in 
contradiction with Eq. (BQ). Therefore, we conclude that 
i?a(A) ~ exp (— AA^/") as A ^ 00. This behavior implies 
that the distribution Qa{Ya) vanishes according to 



]c.{Ya) 



exp 



BYa 



(31) 



as Ya 0. Therefore, the distribution Qa{Ya) has an es- 
sential singularity at the origin, which completes a count- 
able set of algebraic singularities located at Ya = 
with fc = 2, 3, . . . 

We performed Monte Carlo simulations of the stochas- 
tic fragmentation process. Hereinafter, we present sim- 
ulation results for a representative case oi p = 0.4. The 
data corresponds to an average over 5 x 10^^ realizations. 
Figure |l| shows the probability distribution of the second 
moment. The distribution exhibits pronounced singu- 
larities at Y2 = 1/2 and Y2 = 1/3, while the following 
singularities are less visible. One can verify the existence 
of further singularities by differentiating Q2(Y2). Figure 
2 displays the essential singularity at the origin. 




FIG. 1. The distribution of the second moment, (32(^2) 
versus Y2 = ^ . xf, from numerical simulations with p = 0.4. 
The data represents an average over 5 x lO^'^ realizations. 
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FIG. 2. The small size tail of the distribution of the second 
moment. Shown is (32(^^2) versus 1/12- 



IV. EXTREMAL CHARACTERISTICS 

Extremal properties can be viewed as an additional 
probe of sample-to-sample fluctuations. Moreover, they 
are interesting on their own as they arise in many prob- 
lems of mathematics, physics, and computer science 
p2|-^. The largest fragment is an important extremal 
characteristic. Obviously, when x > 1/2, the size dis- 
tribution L{x) of the largest fragment equals the length 
density, L(x) = ^"(2;) = 2qx~^^. In the complementary 
case of a; < 1/2, L{x) satisfies 

1 

L{x) = 2qpL. (^^) + 2p I ^ L (^^) (32) 

l-x 

1-x 

where L-{u) = dv L{v). The first term on the right- 
hand side of Eq. ( p2| ) describes the situation where the 
unit interval is fragmented into two intervals of lengths 
X and 1 ~ x, and where the smaller fragment is stable 
and the larger fragment is unstable (hence the factor 
qp). The latter L_ factor guarantees that subsequent 
fragmentation of the unstable interval does not lead to a 
longer fragment. If one of the first generation fragments 
is shorter than x, then only the longest first generation 
fragment contributes, which leads to the second term on 
the right-hand side of Eq. (^2|). The next term describes 
the situation where both of the first generation fragments 
are longer than x, so the longest fragment can result from 
breaking any of the two fragments. The factor L_ guar- 
antees that the longest fragment of length x comes from 
the corresponding first generation fragment, and the fac- 
tor guarantees that both of the first generation frag- 
ments are unstable. 



Figure 1^ shows that L{x) is discontinuous at x = 1/2. 
This can be understood by noting that L{x) obeys dif- 
ferent equations for x > 1/2 and x < 1/2, and hence it 
loses analyticity at the boundary. Furthermore, we can 
insert the known result, L(x) = 2qx^^P for x > 1/2, into 
the right-hand side of Eq. ( p2| ) and solve for L{x) in the 
interval 1/3 < x < 1/2. Then we can use this solution 
to determine L{x) in the interval 1/4 < x < 1/3, and so 
on. Hence, L{x) should possess an infinite set of singu- 
larities ai X = 1/k, which become weaker as k increases. 
One can also understand why L{x) is discontinuous at 
X — 1/2 by considering the p = case where the distri- 
bution becomes a step function L{x) — 29{x — ^). 




FIG. 3. The size distribution of the longest fragment. 

Consider now the complementary extremal charac- 
teristic - the shortest-segment size distribution S{x). 
Clearly, S{x) = for a; > 1/2. If x < 1/2, one easily finds 
S{x) = 26'(i — x) in the special case of p = 0. To proceed 
in the general case, we first note that, if the unit interval 
is divided into N fragments, the shortest fragment must 
obey Xmin < 1/N . Hence 1/3 < a; < 1/2 implies that the 
unit interval has been divided just once, i.e., both of the 
first generation fragments are stable. This shows that 
S{x) = 2g2 when 1/3 < a; < 1/2. Finally, for a; < 1/3 
the shortest size distribution S{x) obeys 



S{x) = 2q^ + 2qpS- 

l-x 

-h 2p 



y \y 



\-x, 

q+pS^ 



(33) 



where S+{u) = dv S{v). The first term on the right- 
hand side of Eq. ( ^^ describes the situation where both 
of the first generation fragments are stable. The second 
term corresponds to the case where the smaller first gen- 
eration fragment of length x is stable while the longer 
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fragment is unstable, with the 5*+ factor ensuring that 
subsequent fragmentation of this longer fragment does 
not produce a fragment shorter than x. The last term 
describes various situations that are possible if both of 
the first generation fragments are longer than x. 

We cannot obtain an analytical expression for S{x) 
over the entire length range, because in every interval 
(^, a different analytical expression holds. In prin- 
ciple, however, one could determine S{x) recursively. For 
instance, we already know S{x) in the first two regions. 
Inserting those expressions into Eq. ( |33| ) yields 



1-x 



In- 



2x 



(34) 



in the third region 1/4 < x < 1/3. Clearly, S{x) pos- 
sesses an infinite set of singularities at x = 1/fc. 

Figure ^ shows S{x) for x < 1/2 and p — 0.4. One 
can see the plateau region l/3<a;<l/2, and the value 
of S{x) in this region agrees with the theoretical pre- 
diction S{x) — 2q^. Notice that the divergence in the 
small size limit is consistent with the power law behav- 
ior, <5'(a;) ~ x~^ , as a; — !■ 0. Using this observation, we 
insert it into Eq. (|33|), and find a self-consistent value of 
the exponent 6 = 2p. Numerical simulations show that 
S{x) slowly approaches the predicted behavior for the 
case p = 0.4 (see Fig. 5). 

In deriving the relation S = 2p, we have implicitly as- 
sumed that the shortest-segment size distribution is non- 
singular. This is indeed the case whenp < pc- For p > pc, 
however, the distribution S{x) should additionally con- 
tain the singular component 



'S'smg(a;) = AS{x) 



(35) 



with A = 1 — reflecting that with finite probability, 
the total number of fragments is infinite, see Eq. (|2^). 

A more direct way to derive the same result is to 
note that A, the probability that Xmin = 0, satisfies 
A = 2pqA + p'^[l ~ {1 - A)2]. Indeed, the first term 
describes the situation when exactly one first generation 
fragment is unstable while the second term describes the 
situation when both of the first generation fragments are 
unstable. By solving the above equation we find two 
solutions, A = and A = (2p — l)/_p^. The first so- 
lution applies when p < pc] the second solution applies 
when p > Pc and agrees with Eq. (35). In order to in- 
vestigate the small-size asymptotics of S{x), we write 
S{x) = AS{x) + S{x) and assume that the continuous 
part follows the power law behavior, S{x) ~ x^^ as 
X — > 0. Substituting this into Eq. ( p3| ) and balancing 
the dominant terms yields d = 2q. To summarize, differ- 
ent behaviors characterize the small size tail of S{x) 



Of ^ f P< 1/2; 

^^''>^\x-^'i p>l/2. 



(36) 




FIG. 4. The size distribution of the shortest fragment. 



10' 



10' 



10" 



simulation 
slope=-0.8 
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FIG. 5. The small size tail of the size distribution of the 
shortest fragment for p = 0.4. A line of slope —0.8 is shown 
for reference. 



V. APPLICATION TO RANDOM SEQUENTIAL 
ADSORPTION 

Random sequential adsorption (RSA) processes l39| ] 
have been applied to a wide range of chemical, biological, 
and physical processes. Examples include bindingof pro- 
teins to surfaces [^0|j 4l| , g enome sequencing [T5|-[r^, and 
granular compaction ||42| , [43| . In one dimension, fc-mers 
are deposited (with a uniform rate) onto a linear lattice, 
and the deposition events are successful only when all k 
sites are empty. Eventually, the system reaches a jam- 
ming state where no further deposition events are possi- 
ble. Basic quantities of interest are the jamming density. 
Poo, and its dependence on the initial concentration po, 
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as well as the gap size distribution. 

Adsorption can be viewed as a fragmentation process, 
and the above stochastic fragmentation process general- 
izes RSA to situations where the gaps between the ad- 
sorbed particles may become passive with probability p 
after each deposition event. Stable gaps can no longer be 
deposited onto, or in other words, fragmented. In the fol- 
lowing, we study this RSA problem both in discrete and 
continuous space. We show that the limiting density of 
passive gaps is proportional to the fragment length den- 
sity obtained above. In the interesting limit of vanishing 
initial concentrations, po ~^ Oj the final jamming density 
vanishes according to p^o ~ q < 1/2, and as poc ~ Po 
if q > 1/2. Hence, the jamming density is significantly 
enhanced over the initial density only if p exceeds the 
critical value pc = 1/2. 



A. Discrete Space 

Let us consider random sequential adsorption in one 
spatial dimension were each deposition event creates two 
new smaller gaps. We assume that each of these gaps 
remains active with probability p, while with probability 
q = 1 — p it becomes passive, i.e., adsorption events no 
longer occur on this gap. Initially, the system consists 
of randomly distributed monomers with density po, and 
all gaps are active. Then, r-mers are deposited with a 
temporally constant and spatially homogeneous rate, set 
to unity without loss of generality. The deposition events 
are successful if and only if all r sites are empty. Below, 
we consider the monomer case (r — 1). 

Let Ak{t) {Pk{t)) denotes the probability of finding an 
active (passive) gap of size k at time t. These gap proba- 
bility distributions satisfy the following master equations 



dAkjt) 
dt 

dPkjt) 
dt 



2p J2 Anit)-kAkit), 



2q J2 

rn— 



(37) 



The loss term reflects that deposition in each of the 
empty sites destroys the gap, and the gain term reflects 
the fact that two smaller gaps are created in each de- 
position event. The prefactors of the gain term equal 
the corresponding production rates. It is easy to ver- 
ify that the evolution equations conserve the total length 

Efc(fc + i)(Afc + Pfc) = i- 

Let us consider systems that initially consist of ran- 
domly distributed monomers with all gaps being active, 
i.e., Ak{0) ^ pI{1~ Po)'' and Pfc(O) = 0. Since Pk is 
enslaved to Ak-, we derive the latter quantity first. The 
linear loss rate suggests the exponential ansatz 



Ak = aP'' 



(38) 



with the initial values a(0) — pg and /3(0) — 1 — po- 
Substituting (|38|) into the rate equation yields 



df3 d 

—— = —L>, — in Of 

dt dt 



2p 



1-13 



(39) 



The first equation yields P{t) = (1 — po)e~*, and the sec- 
ond equation can be conveniently solved by changing the 
time variable from t to P using (3dt = —dp. This trans- 
forms the second equation into ^Ina = —2p{l — P)~^ . 
Integrating this equation subject to the above initial con- 
dition yields a = pl''{l — PY^ ■ The time-dependent gap 
distribution is therefore 



Ai 



Pf^p\ /3=(l-po)e- 



(40) 



Next, we study the final jamming density which can 
be obtained from the active gap distribution by inte- 
gration of the overall deposition rate over time, i.e., 
Poo-Po = C Efc kAk{t) - dtaP{l-P)-\ Again, 
it is useful to transform t to p. Evaluating the integral 
gives 



Poo 



(p^'~2gpo)/(l 
Po ln(l/po) 



2q) 



9^1/2; 
9=1/2. 



(41) 



It is easy to verify that poo — > 1 in the limit g — > 0. In the 
more interesting limit when the system is initially almost 
empty, i.e., when po — > 0, the following leading behaviors 
emerge 



Po' 9 < 1/2; 

Poln(l/po) (7=1/2; 
Po q> 1/2. 



(42) 



In other words, the final coverage depends algebraically 
on the initial coverage in the range q < 1/2. In this case, 
the adsorption process can be viewed as "effective" since 
the increase in density is significant. Otherwise, the fi- 
nal density is proportional to the initial density po. The 
critical point is marked by a weak logarithmic increase 
in the jamming density. 

We turn now to the distribution of passive gaps which 
can be obtained by integrating the master equation (|3^ ) 
with Ak{t) given by Eq. (^0|). In order to compare the 
gap-size distribution with the fragmentation case, we fo- 
cus on the limiting {t oo) distribution of passive gaps, 
which reads 



Pfe(oo) = 2qpo'^ 



l-po 



dp (1 - P)^p~^p'' 



(43) 



0, 



In the limit of almost empty initial conditions, po — 
where the average gap size diverges, we find (fc) = 

(1 - Poo)/Poo ~ P^^ 

Hence, the most interesting behavior emerges in the 
scaling region po — *■ and fc — > oo with the scaling vari- 
able ^ — kpo kept fixed. In this region the limiting gap 
distribution ([43|) can be rewritten in the scaling form 



P,(^) = 2gp2«A;-2pr(2p,0 



(44) 
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where r(a. ^) = f^dxx'^ ^ denotes the incomplete 
gamma function. Equation (jij) shows that the gap size 
distribution behaves algebraically as long as the size of 
the gap does not exceed the average initial size k* = Pq^ , 
while for larger gaps the size distribution is suppressed 
exponentially 



1 2qpok-'^e-^/^' k :^ k* 



(45) 



These two expressions are indeed of the same order, 
P ~ Pq, in the vicinity of the crossover point k ^ k* . 
In general, we find the algebraic behavior Pfc(oo) ~ fc^^^ 
in the limit po — > 0, in agreement with the stable frag- 
ment length density of Eq. (^) . This agrees with intuition 
since in the limit po ^ the stochastic RSA is a discrete 
counterpart of the stochastic fragmentation. 

The above treatment can be generalized to the dimer 
(r — 2) case and even to the general r-mer case. Al- 
though these solutions become very cumbersome as r 
grows, the asymptotic behavior found for the monomer 
case including the scaling form of Pk{oo) and the jam- 
ming density p^o are not altered. 



B. Continuous Space 

The continuum limit where particles of unit length are 
deposited irreversibly onto a line can be obtained from 
the discrete r-mer case by taking the limit r oo and by 
redefining the time variable rt t and the initial density 
rpo A. Then, the densities of active and passive gaps 
of size X evolve according to 



dA{x,t) 

dt 
dP{x,t) 

dt 



2p 
2q 



dyA{y,t) - 0{x - l)(a;- l)A{x,t), 



dyA{y,t), 



(46) 



x+l 



where 9{x) is the step function. The initial conditions 
read A{x,0) = A^e"^^ and P(x,0) = 0. 

We first derive the distribution of active gaps of lengths 
a; > 1. Equation ( p6| ) suggests that it remains exponen- 
tial throughout the evolution, i.e.. 



A{x, t) = $(i) exp [-Aa; - (x - l)t\ 



(47) 



Substituting this exponential form into Eq. (^) yields 
^ln<i>(t) = 2pe~^~''' /{\ + t). Integrating this differen- 
tial equation subject to the initial conditions ^(O) = A^ 
gives the time dependent prefactor 



$(t) = X^\\ + t) 



2p 



exp 



-2p 



dr. 



The jamming density can be obtained by 
grating the total deposition rate over time, pr 
Jp dt dx {x — l)P{x,t), which yields 



(48) 



inte- 



Po 



dt 
t^ 



exp 



-2p I dr- 



(49) 



When p = 1, this expression agrees with the jamming 
density of the parking model with and without disorder 
p^ , p5[ . Independent of the probability p, the approach 
to the jamming state follows the classical t~^ law [p9[ 



Po 



p(t) ~ A2«e-^(A + <)-i 



(50) 



Additionally, the leading behavior in the limit of dilute 
initial conditions (A — > 0) can be evaluated and the be- 
havior found in Eq. (|4^) generalizes to the continuum 
limit 



A29 g<l/2; 
Aln(l/A) q=l/2- 
X q> 1/2. 



(51) 



The limiting passive gap distribution is found by inte- 
grating the rate equation ( p6| ) using the active gap dis- 
tribution 



Pooix) = 2q 



dt 

X + t 



$(t)e 



-A-(A+t)a: 



(52) 



with ^(t) given by Eq. (^). In the limit of dilute initial 
conditions (A 0) one can simplify the integral on the 
right-hand side of Eq. (^2|) to find the following extremal 
behaviors of the gap distribution 



( 2qXx-^e-^'' 
,{x) ~ <^ 2gr(2p)A29x-2p 

[ 2ge-2p7EA29 In(^) 



X > A~^ 

1 < a; < A"^ 

a: < 1: 



(53) 



where 7^ = 0.577215 . . . denotes the Euler constant. The 
first two asymptotics in Eq. ( |53| ) are straightforward ex- 
tensions of the corresponding behaviors in the lattice 
case, and the last line of Eq. ( [53| ) has been derived from 
Eq. (H) using dT{l~e-^)/T = lnT + -fE + 0{e-^). 

Finally, we note that even in the long-time limit the 
density of active gaps does not vanish for sufficiently 
short gaps, a; < 1. One can determine A^(x), and more 
generally A{x,t), by employing an elementary relation 
between the densities of active and passive gaps, namely 



A{x,t) = A{x,0) + -P{x,t). 



q 



(54) 



This relation immediately follows from the master equa- 
tions (46), and it clearly holds for arbitrary t as long as 
X < 1. By combining (^3|) and ( |5^ we find 

A^{x) ~ 2pe-2P7^ X^'^ In j , (55) 

which applies if A ^ 1 and a: ^ 1. 

Therefore, stochastic fragmentation processes can be 
naturally extended to adsorption processes, and apart 
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from numeric prefactors the algebraic fragment-size dis- 
tribution is reproduced in the hmit of empty initial con- 
ditions. The phase transition underlying the branching 
process has an interesting implication. The jamming den- 
sity is significantly larger than the initial density only 
when p > 1/2. The super-critical nature of the under- 
lying branching process allows for an infinite number of 
fragments produced from a single fragment, and this ex- 
plains the enhanced jamming density in the stochastic 
RSA process. 

Although the fragmentation and the adsorption results 
are closely related, we have used two complementary ap- 
proaches to obtain them. In the former case, it was con- 
venient to bypass the distribution of unstable fragments 
and solve directly for the final stable fragment distribu- 
tion, while in the latter case, it was more natural to study 
the entire time dependent behavior of both distributions. 
A more complete treatment of the fragmentation process 
is of course possible using a continuous time formulation 
which leads to rate equations similar to Eq. (Eq). 



VI. SUMMARY 

We have studied a class of stochastic fragmentation 
processes, where fragments may become stable ( "frozen" ) 
after each fragmentation event. We have found that in 
general, these processes are characterized by an algebraic 
small-size divergence of the fragment size distribution. 
This behavior is robust as it holds for size dependent 
fragmentation densities and fragmentation probabilities, 
as well as in dual adsorption processes in both continuous 
and discrete space. The corresponding power-law expo- 
nents can be tuned by varying the fragmentation proba- 
bility, and the entire range allowed by mass conservation 
may be realized. 

While the size density can be determined analytically, 
additional statistical measures of fluctuations are more 
difficult to handle. Nevertheless, we have shown that mo- 
ments of the distribution exhibit large sample-to-sample 
fluctuations, and hence, knowledge of the entire distribu- 
tion of observables is needed to characterize the system. 
Additionally, the distribution of the moments and of ex- 
tremal characteristics, such as the longest and the short- 
est fragments, possesses an infinite set of singularities. 

Lack of self-averaging is important in practical applica- 
tions such as utilization of DNA segmentation for com- 
parison of genomes of different species. Great care is 
clearly needed in comparative analysis of the segment 
length distributions as observed deviations between seg- 
ments may be actually statistical rather than biological. 

We thank Michael Zhang for valuable discussions. This 
research was supported by DOE (W-7405-ENG-36), NSF 
(DMR9978902), ARO (DAAD19-99-1-0173), NIH, and 
the CSHL Association. 
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